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ABSTRACT 


Damage reduces the flexural stiffness of a structure, thereby altering its dynamic 
response. Considerable effort has been put into obtaining a correlation between the changes in 
modal parameters and the location and amount of the damage within the structure. Most 
numerical research employed elements with reduced beam stiffness to simulate damage in the 
beam. This approach to damage simulation neglects the non-linear effect that a crack has on the 
structural dynamic response. In the present study, finite element modeling techniques are 
utilized to directly represent an embedded crack. The results of | the dynamic analysis of the 
present model are then compared to the results of the dynamic analysis of the reduced modulus 
finite element model. Different modal parameters are investigated to determine the most 
sensitive indicator of damage and its location. Nonlinear effects, such as crack closure and 
opening, of an embedded crack on the structural dynamic response were also studied from 
transient nonlinear analysis. The modeling technique is then applied to sandwich composite 
beams with simulated delamination to investigate damage detection techniques through the use 


of damping caused by frictional dissipation of energy on the crack surface. 
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I. INTRODUCTION 


The ability to monitor a structure and detect damage at 
the earliest possible stage is a very important asset in all 
realms of civil, mechanical, and aerospace engineering 
disciplines. Many existing methods of damage detection use 
various non-destructive techniques such as visual, acoustic 
emission, ultrasounds, magnetic fields, radiographs, eddy- 
currents, C-scan, and vibration thermographs. Almost all of 
these meenode are inherently localized and require prior 
knowledge of the possible damage location in order to avoid 
wasting time during useless data collection. Using these 
methods for large scale or global areas can be sae: 
expensive as well as require large amounts of labor and 
time. Therefore, an alternate approach that incorporates the 
evaluation of a single global parameter or combination of 
global parameters to detect, localize, and monitor damage 
would greatly enhance the sAeaes and efficient maintenance 


of structures. 


More recent work in the area of global damage detection 
techniques has involved the use of relative shifts or 
changes in the modal parameters. Modal parameters such as 


modal frequency, modal damping, and mode shapes are a 











function of the mass-inertia and elastic iittala of the 
seeaeHiEe: Therefore, the dynamic response of a system can 
be a sensitive indicator of change in the integrity of the 
system’s elastic structure. Damage in any form reduces the 
local flexural stiffness in the Vicinity of the damage. 
Reduced suitensss leads to a decrease in modal Srenaeneiee 
and an increase modal damping sospnictenrs as well as 
changes in the corresponding mode shapes. 

Damage can affect the dynamic response of a structure 
either linearly or nonlinearly. Changes in the dynamic 
response of the structure subject to linear damage can be 
related to uniform changes a the geometry or the material 
properties of the structure. However, in most cases damage 
is predominantly in the form of micro-scale or macro-scale 
cracks or voids. Such cracks or voids behave nonlinearly 
under dynamic excitation as aia the case where embedded 
cracks subsequently open and close. In Some nin nese cases, 
the global effect of damage on the dynamic response of the 
structure can be considered to be small and local in nature. 
In such cases, the damage can be modeled linearly. However, — 
in the case siere. tie nonlinear local behavior significantly 
affects the global structural response, the damage can no- 


longer be modeled linearly. 








Most of the previous research nee concentrated on 
linearized modeling of damage. In these linear eigenvalue 
analyses using the finite element method, using either 
reduced material anspentiee or reduced geometric dimensions 
within the damaged section simulated damage. The, focus of 
these studies was aimed at establishing a modal parameter or 
method of manipulating modal parameters to obtain the 
highest sensitivity to damage in the structure. However, in 
linearizing crack damage in this fashion most researchers 
must reduce the flexural stiffness in the damaged elements 
in large amounts to produce relatively small changes in the 
modal parameters. Little thought has been centered on the 
effects of linearization of local crack response or in the 
direct representation of an embedded local crack in a finite 
element model. 

One of the objectives of this study is to evaluate the 
previous modeling technique, which smeared a crack into a 
reduced material or geometric properties. To this end, the 
present study includes an embedded crack modeled directly 
instead of smearing the crack through reduced material or 
geometric properties. Linear eigenvalue analysis was 
performed for this model to compare the present modeling 


technique with previous smearing techniques used in damage 








detection research. In this study, the crack behavior was 
linearized. 

“The second objective was to investigate the modal: 
parameter or combination of modal parameters that was more 
sensitive to a local embedded crack using the present 
modeling technique. For example, modal frequency, mode shape 
displacement, mode shape curvature, etc. were examined for a 
possible modal parameter, which would provide accurate 
information regarding the presence of a crack and its 
location. 

The final objective was to sane the local nonlinear 
behavior (such as crack opening and closure) of a crack 
using transient nonlinear analysis. The effort was made to 
evaluate the linear eigenvalue analysis conducted in the 
first portion of the study. Furthermore, the transient 
nonlinear analysis was utilized to investigate damage 
detection and localization through a new modal damping 


(energy dissipation) parameter. 











II. LITERATURE REVIEW 


There is quite a large amount of literature generated 
on global damage detection and localization by use of 
changes in modal parameters. The majority of literature 
pouaewes for the current research has been focused on damage 
modeling techniques, the moaaa parameters being utilized for 
damage detection and localization, and the level of 
sensitivity to the presence of damage. 

The amount of literature related to damage detection 
using relative shifts in modal frequencies is quite large. 
Most results summarize a relatively low sensitivity of 
frequency shifts to damage and the need for very precise 
measurements for damage detection. This is mainly due to the 
reduction in flexural stiffness being local in nature, 
thereby resulting in modal frequency shifts that are hardly 
detectable with the current test equipment. Spyrakos (1990) 
found a modal ppeaueney shift of less than 5% corresponding 
to a 50% reduction in flange cross-sectional area. Rehm 
(1987) also found large-scale reduction in eeseeeeeeeieaal 
area (over 30%) only produced small and less detectable 
shifts in the modal frequency (less than 1%). Because a. 


shift in modal frequency is a global parameter of the 








structure, it provides no clear and precise means of 
identifying more than damage existence. Damage of the same 
size but in different locations of the structure may produce 
the same amount of frequency change. This has led many 
researchers to investigate the use of more sensitive modal 
parameters to detect and locate damage. 

Formulas.and algorithms wenng sieeiaceneae mode shapes 
and curvature mode shapes have been used to better locate 
the damage. Yuen (1985) in his paper investigated the first 
five mode shapes of a damaged cantilever beam, varying the 
location and the degree of damage within the beam. He 
modeled the crack as a reduction in the modulus. of 
elasticity. He borrowed this procedure of reducing the 
modulus of elasticity from Adams (1978), who reduced the 
aoaae of elasticity to zero in the damaged area. Yuen 
(1985) found that a reduction in the modulus of elasticity 
of 50% produced a small but measurable change in both the 
first modal frequency and mode shape. He further added that 
the expected change in vibration amplitude was measurable 
and would demand the same order of accuracy as the frequency 
measurement. In this study, he failed to mention the exact 
order of accuracy required. In a study conducted by Pandey 
(1994), changes in mode shape ‘eu were used to 


identify and locate damage in a simply supported beam and a 








cantilever beam model. Pandey (1994) reverted to the common 
practice of reducing the modulus of elasticity to model the 
damaged section. However, he showed that damages detected by 
changes in mode shape curvature were localized to the region 
of damage while changes in the mode shape displacement were 
not localized. 

The carrying precedence in the majority of articles 
reviewed was crack damage continuously being modeled as a 
eduction in the flexural stiffness (EI) of the damaged 
elements. In almost all cases, there was no considerable 
effort provided to explain the linearization of dynamic 
response for neither crack damage nor sufficient evidence to 
support the amount of material property or geometric 
property reduction used in the analysis. Following this 
procedure, there is no clear relationship between the 
geometrical dimensions of the crack and the amount of 
reduced flexural stiffness incurred within the structure. 
This method also provides no exact correlation between the 
damaged section’s effect on the global and local dynamic 
response of the structure among the different modes of 
vibration. Therefore, there is a need to clarify the 
differences between the linearization in dynamic response of 
embedded local crack damage and the direct representation of 


embedded local crack damage in a finite element model. 
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IIL. FINITE ELEMENT MODELING AND ANALYSIS 


A finite element model was generated using beam 
elements with displacement degrees of freedom only. The 
detailed ee ation is given by Kwon and Bang (1997) and is 
therefore omitted in the following sections. As seen in 
Fig. 1, each element has six degrees of freedom. There are 
four axial displacement dof’s at the corner points ane two 
transverse displacement dof’s at the endpoints of the 


element. 





a5” | aio x 
Figure 1. Four-Noded Beam Element with Six Dof's 


The displacement field of the element is 


wf”) 











where [N] is the linear shape function matrix and {d*°} is 
the vector of nodal displacements. The axial normal strain 


can be written as 


ON N 
a inc — Hu; a Pier 5: ee (2) 
ax ox Ox Ox Ox 


and the shear strain is 


The element stiffness matrix was obtained by minimizing the. 
total strain energy which contains both bending and 


transverse shear energy. 
[Rol = LR diet TR Ie (4) 


The subscripts ‘b’ and ‘s’ indicate bending and transverse 
shear, respectively. The bending and transverse shear 


stiffness matrices are 


T 


(B,} E4B, jdxdy | (5) 


Ks ]= 
Ik: ]= 


iB, , E{B Saxdy (6) 


10 








where E and G are the elastic and shear modulii of the beam 
and the vectors {Bp} and {B,;} are derived from the strain- © 


displacement relationship given below. 


c, = {B,}a°} - (7) 
vy = 1B, }a°} (8) 


By carrying out the intégration, the bending and transverse 


shear stiffness matrices are obtained. 


2 1 0-2 -1 0 
i 2 @ <1 =2 0 

1] Eni 0 0 0 0 0 0 

KS ]=F jako 2 1 0 
1-20 1 2 0 
0 000 0 0 

and 


Gi? -GP? 2Glh Gl’ -Gl’? -2Glh 

-G? GI -2Glh_ -GI’— GI’ 2Ghh 

.) 1 | 2Glhk -2Glh 4Gh? 2Glh -2Glh -4Gh’ 
IKS Ah 2 2 2 2 (10) 

*“ 4lh| GP? -Gl 2Glhh_ = Gl -Gl?  —2Glh 

-G? Gl’? -2Glh -Gl’? Gl’ — 2Gh 

~2Glh 2Glh -4Gh? -2Glh 2Glh 4Gh’ 


tae 








The element stiffness matrix, which is obtained by adding 
the bending and transverse shear stiffness matrices, can be 


expressed in the following form 


(a, +2a,) (-a,+a;) a, (a-2a) (-a,-a;) -a, 
(-a,+a;) (a,+2a;) -a, (-a,-a;) (a,-2a,) a, 


i ]-| 


(-a,—a,) (a,-—2a,) -a, (-a,+a;) (a,+2a;) a, 


(11) 


a,—2a,) (-a,-a;) a, (a,+2a,) (-a,+a,;) -a, 


where each symbol denotes 


Gl Gh Eh G 
ee ee ac aaa a, 61 ary 


In the mass matrix, a lumped mass was used where the mass 


was concentrated at the six nodal degrees of freedom. 


] 
0 
phib| 0 
4:10 
0 
0 


IM‘ |= (12) 


S20O 0 Oro 
2Z2OONOCO 
o2Oor CoO CO Oo 
OF OOO CO) 
NCO O CO Oo 


The phlb term is the element mass and is distributed among 


the four axial displacement dof’s and two transverse 


displacement dof’s. 


i2 





Once the system mass and stiffness matrices are 
computed for the beam, they are substituted into the 


equation of motion of an undamped vibrating system. 
[M]ik}+ [Kix}=0 (13) 


After applying the proper boundary conditions, the above 
equation is ‘solved using the eigenvalue equation given 


below. 
(K-27 IMJjo, =o : a4) 


The term Az; represents the modal frequencies and Qj; 
represents the vector of corresponding displacement mode 


shapes. 


Although damage can be detected from the change in 
displacement mode shapes, finding the location of the damage 
is very difficult. Pandey (1994) found that changes in 
displacement mode shapes are not localized to damage areas. 
Therefore, other parameters which are more sensitive to the 
location of the damage must be utilized. Mode shape 
curvature is related to the flexural stiffness of a beam 


cross-section by the following relationship: 


or-M (15) 


L3 








where ®” is the curvature at the section, M is the Bending 
moment, and EI is the flexural stiffness. From the 
displacement mode shapes, curvature mode shapes were 
obtained numerically by using a weighted central difference 


approximation. 


D ia (x, ~ Hi )- ® (x4 ~ Xj. )+ OD 5, (Ki ~ x,) 


(16) 





Pandey (1994) successfully showed the benefits of using mode 
shape curvature to detect and locate damage in simply 
Supported and cantilever beam, where damage was simulated by 


a 50% reduction in the modulus of elasticity. 


14 





IV. FINITE ELEMENT MODEL COMPARISON 


A. MODEL 

The material pronerries for an aluminum alloy (1100- 
H14) were used to generate a mesh for a 48cm x 2.54cm x 
1.27cm beam. The beam was modeled as a single row of ake 
elements with the majority of elements concentrated around 
the damaged section. A 2.4 cm, through the width, axial 
crack located at the midpoint of the beam was Simulated by 
reducing the modulus of elasticity in the elements of the 
damaged section. Another beam was generated, having two 
rows of eighty elements with the majority of elements 
concentrated around the damaged section. A total of 160 
eoaced elements were used for this model. Transverse 
displacement degrees of freedom ae coupled along the beam 
axis. Axial displacement continuity was enforced along the 
interface of the top and bottom row of elements. An 
embedded, through the thickness, crack was simulated by 
uncoupling the axial and transverse degrees of freedom 
between the top and bottom row of elements. The finite 


element models are included in Figs. 2 and 3. 


15 
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B. RESULTS 


The Sieeng sue problem was solved for the finite 
element model of the undamaged beam and the beam with an 
aipedaed: eae. The first five non-zero modal frequencies 
and corresponding mode shapes were obtained for the damaged 
and undamaged beam using the embedded crack wieder, Only 
transverse displacement dof’s were considered for the mode 
shapes. The absolute shift in the modal frequencies due to 
the presence of the crack was calculated and listed in Table 
1. Once the modal frequency shifts were suena aso: the 
aaniuG, 62 elasticity in the damaged elements of the other 
finite element model were adjusted until the corresponding 
modal frequency shifts were approximately equal to those 
experienced by the finite element model with the embedded 
crack. Simultaneously matching the first five non-zero modal 
frequency shifts was impossible. However, by reducing the 
modulus a very small amount (7x10°'3), the first non-zero 
modal frequency shift was matched closely as seen in Table 
Bis As seen in Table 3, reducing the modulus’ a very large 
amount (75%) matched the second and the fourth non-zero 
modal frequency shifts between the embedded crack model and 


the reduced modulus model. 
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Table 1. Frequencies for Embedded Crack Model 






[a _[ 28-6688 | 28.6088 | 0 
a Te ezat T8652 | a5 
[>a [250.6039 [249.1599 | 1.4440 
(s|369.9797 [369.9429 | 0368 


Undamaged amaged [A-A! | (Hz) 
AX (Hz) : (Hz) 








Table 2. Frequencies for 7e-7% Reduced Modulus Model 











Undamaged | Damaged | |\-A'| (Hz) 
(Hz) dX! (Hz) 


[2 [78.2601 | 78.2601 | ~0 
[4 [247.7349 | 247.7349 | =O 
[5 | 364.2232 | 364.2231 | 0001 





Table 3. Frequencies for 25% Reduced Modulus Model 











(Hz) XN’ (Hz) 
[2 [78.2601 | 78.1053 | 0.1548 
[4 [- 247. 7349_| 246.2816 | 1.4534 
[s_[ 364.2232 [336.8181 | 27.4050 


An average of the first five modal 


frequencies 


Was 


calculated for the embedded crack model and corresponded to 


the average frequency change of a model with a 6% 
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decrease 














in the modulus of elasticity in the damaged section. The 
values for the modal frequencies are listed below in Table 
4. Although the average change in the first five non-zero 
modal frequencies matched the average change of the embedded 
crack meee the individual changes in modal frequency nor 
their corresponding changes in displacement mode shapes 


matched the changes experienced by the ‘embedded crack model. 


Table 4. Frequencies for 94% reduced Modulus Model 


Undamaged A (Hz) | Damaged A’ (Hz) [A-A’ | (Bz) 
28.6075 28.4938 ey, 
a 78.2601 78.2568 .0033 
rT 3Cid 151.6761 151.2246 4515 
ae 








Although the modal frequency shifts caused by the crack 
presence could be matched sees cay, matching isis 
corresponding changes in the displacement mode shapes was 
rather difficult. Even numbered displacement mode shapes 
matched with close accuracy when matching their 
corresponding change in modal frequency. However, the 
displacement mode shape changes for the first mode shape and 
all other odd numbered modes shapes could not be matched to 
the changes experienced by the embedded crack model by any 


alteration of the modulus of elasticity in the damaged 


section. Altering the modulus of elasticity in the damaged 
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section had a more site effect on the change in the odd 
numbered displacement mode shapes than the embedded crack 
“model for the given damage scenario. The magnitude of 
displacement change in the reduced section modulus model was 
also found to be significantly lower than the displacement 
change found in the embedded crack model when the first 
modal frequency. change was matched. It was very clear that 
ere existed some difference between the models that could 
not be corrected by any puseievebie alteration in the 


modulus of elasticity in the damaged section. 
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V. LINEAR DYNAMIC ANALYSIS 


A. MODEL 

The same material properties for the pen alloy 
(1100-H14) beam used in the comparison of reduced material 
property damage and the embedded crack damage models were 
used to generate the mesh for the first study in the Linear 
Dynamic Analysis. Having acknowledged the differences 
between the dynamic response of an embedded crack model aa 
a reduced modulus of elasticity model, the modal 
frequencies, mode shape displacements, and mode_ shape 
curvatures were calculated for a free-free beam with and 
without a 2.4 cm embedded crack wen the embedded crack 
model. Each modal parameter was compared for sence vies tO 


damage detection and localization. 


The analysis was further carried out on a finite 
element mesh representation of a sandwiched composite beam 
with the following dimensions 48 cm x 4 cm x .8 cm. In the 
finite element model, Syntac 350C foam of .6 cm cna eknese 
(Young’s Modulus=2.21 x 10° N/cm? Poisson's Ratio=.35, 
Density=9.61 x Lo kg/cm, Modulus of Rigidity= 8.18 x 10” 
N/cm?) eae reinforced by two .1 cm thick layers of Glass 


Reinforced Plastic (Young’s Modulus=2.07 x 10° N/cm* 
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Poisson’s Ratio=.342, Density=3.06 x alk 6 te kg/cm’, Modulus Of 
Rigidity= 8.97 x 10° N/cm*). The finite element model of the 
composite beam is included in Fig. 4. The laminae were 
modeled di een eeceny joined at the interface except in the 
damaged section, where delamination was simulated. The size 
of the crack was varied between 2.4 cm and 3.6 cm. The 
location of the crack was also varied between the beam 
midpoint and 10 cm away from the midpoint. Three separate 
boundary conditions were applied to the model. (Free-Free, 


Simply Supported, and Cantilever) 









Patni Crack 
TS 
PEE, CECH ie FOAM 






Figure 4. Finite Element. Mesh for Composite Beam with Crack 


B. RESULTS 

The data collected from the embedded crack model was 
used to detect and locate the embedded, through the 
einess, crack. The maximum values of percent relative 
change for each curvature mode shape and displacement mode 
shape as well as the percent relative change in modal. 
frequency are listed below in Table 5. As previously 
mentioned, the modal frequencies are almost insensitive to 


the embedded crack damage. It was also hard to derive the 
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location of the crack from the modal frequency shift due the 
possibility of different cracks producing the same shift in 
modal frequency. Therefore, the analysis was more focused on 
the mode shape displacement and mode shape curvature values. 
In all modes, the maximum values of relative change were 
higher for curvature. mode shapes, displaying a higher 


sensitivity to the presence of damage. 


Table 5. Maximum Percent Change 


[Mode [Percent AX | Maximum percent Ad| Max. percent Ag!’ | 
7 


Plots of the displacement and curvature mode shapes are 















included in Figs. 5 through 14. As seen in the figures, 
changes in curvature mode shapes were more localized to the 


region of damage than changes in displacement mode shapes. 
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First Mode (Omega=28.6688Hz) 
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Figure 5. First Modal Displacement of Embedded Crack Model 
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Figure 6. 





Second Mode (Omega=78.6241Hz) 
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Second Modal Displacement of Embedded Crack Model 
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Modal Amplitude 


Third Mode (Omega=152.8576Hz) 
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Figure 7. Third Modal Displacement of Embedded Crack Model 
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Fourth Mode (Omega=250.6039Hz) 


Modal Amplitude 
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Figure 8. Fourth Modal Displacement of Embedded Crack Model 


27 











Modal Amplitude 


Fifth Mode (Omega=369.9797Hz) 
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Figure 9. Fifth Modal Displacement of Embedded Crack Model 
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Figure 10. First Curvature Mode of Embedded Crack Model 
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Figure 11. 
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Second Curvature Mode of Embedded Crack Model 
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Third Mode (Omega=152.85/6HZz) 
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Figure 12. Third Curvature Mode of Embedded Crack Model 
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Figure 13. Fourth Curvature Mode of Embedded Crack Model 
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Figure 14. Fifth Curvature Mode of Embedded Crack Model 
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Modal response data was collected for three different 
boundary conditions with the size and location of the crack 
varied. The mode shape changes and mode shape curvature 
changes were evaluated for each case with hopes of 
solidifying the change in modal peepones ah terms of crack 
size, location, and boundary condition. 

It was found that the boundary condition applied had a 
less significant effect on the modal response change 
generated by an embedded crack of a fixed size and location. 
This can = a eaees to the fact that the odd numbered 
modes of vibration for a cantilever beam are similar to 
those of a free/free beam due to symmetry and the response 
of a simply supported beam is very similar to the response 
of a free/free beam. 

Increasing the crack size increased the amount of 
change in the modal response, whether considering frequency, 
mode shape displacement, or mode shape curvature. However, 
there was not a direct proportional relationship that could 
be established between the crack size and the amount of 
change experienced. The change experienced by each mode was 
Significantly different for each case of damage. 

The location of the crack in a given mode of vibration 
affected the amount of crack opening and closure (nonlinear 


effect) experienced in the response which relates to the 
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detectability of damage from changes in the modal response. 
From varying the location of the crack, it was found that 
the local nature of the crack’s effect on modal response 
became more significant in modes where the reduced stiffness 
was more evident through crack opening. These locations were 
around the anti-nodal locations of the mode of vibration. A 
plot of the more sensitive and less sensitive crack 


locations is provided below in Figure lo. 


Mode 2 


0.8 | ‘less effect 


0.2 more effect 


S 


Curvature or Amplitude 


Oo 9 
on 
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Figure 15. Modal Sensitivity to Damage Plot 
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VI. NONLINEAR DYNAMIC ANALYSIS 


A. INTRODUCTION 
The preceding analysis has concentrated on the 
linearized effects of crack damage through the use of the 


eigenvalue problem. The Finite Element mesh generated thus 





far did not account for contact between the interior 
surfaces of the crack during a dynamic response. The data 
collected therefore considers the problem of modal parameter 
changes due to the effects of crack damage in terms of a 
stiffness reduction in the location of the crack only. 
Therefore, an analysis of the transient response of a 
damaged beam with consideration for the contact surfaces is 
studied for capturing both the effects of crack surfaces 
impacting and the nonlinear opening and closure of the crack 


on the change in dynamic response. 


B. ANALYSIS AND RESULTS 

To this end, LS-DYNA software was incorporated for the 
analysis of the nonlinear transient Pesponse: A mesh with 
elements of equal length was generated for the same 


dimension composite beam used in the Linear Analysis portion 
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of the current research. Cantilever boundary conditions were 


applied to the beam which is depicted below in Figure 16. 
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Figure 16. Transient Analysis Model 


The location and size of the crack, respectively labeled ‘a’ 
and ‘b’ in the figure above, were eer for the analysis. 
Natural Frequency values for the first mode of vibration 
were very close to the values obtained in the linearized 
analyis. Frictional impact surfaces were created between the 
surfaces of the crack. The friction coefficient was varied 
to evaluate the amount of frictional dissipation of energy 
involved with crack vibration during a transient response. 
Since the exact value for the friction coefficient for a 
crack surface was unknown, values were chosen in an effort 
to relate the lack of smoothness that would be found on the 
crack surface. The embedded crack produces an increase in 
damping which can be captured by a decrease in sucessive 
amplitudes during the transient response of the beam. The 


following logarithmic decrement equation was used to 
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calculate the damping ratio values for two consecutive 


cycles in tip deflection of the beam, 


(17) 





where ¢ is the damping ratio and x, and x2 are successive 
amplitudes. The results From varying the £EVCLLON 


coefficient are included below in Table 6. 

Table 6. Damping Ratio According to Friction Coefficient 
Friction Damping ratio 
Coefficient (1) 


Note: a=24 cm, b=2.4 cm 





The data shows a nonlinear relationship between the friction 
coefficient and the amount of damping present in the 
response of a eonposiee beam with a 2.4 cm crack embedded at 
the midpoint of the beam length. The results from increasing 


the crack size are included in Table 7. 


Table 7. Damping Ratio According to Crack Size 


Crack Size ‘b’ | Damping Ratio 
(cm) 


0.0222 
pO 0281 


Note: w=1,a=24 cm 
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As observed in the preceding table, there was a 21% increase 
in modal damping caused by a 40% increase crack size. The 
effect of crack location on modal damping is summarized in 


Table 8. 


Table 8. Damping Ratio According to Crack Location 


Crack Location ‘a’ Damping Ratio 
(cm) 


14 0.0153 
24 Oc0Z22 
.34 0.0585 


Note: p=1, b=2.4 cm 












Important to the analysis was the nonlinear opening and 


closing of the embedded crack. By plotting the relative’ 


axial and transverse displacements of two nodes in direct 
contact at the midpoint of the crack surface, a clearer 
picture can be drawn about the periods of contact between 
the crack surfaces and the amount of nonlinear dissipation 
of energy through crack surface friction involved when the 
crack surfaces are in contact. The axial and transverse 
displacement plots of two nodes at the top and bottom of the 
crack surfaces (see nodes A and B in Figure 10) during the 
transient response of a composite beam with a’ 4.0 cm crack 
embedded crack at the midpoint of the beam length are 


provided in Figures 17 and 18. 
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x 10 Crack Opening vs Time 


Crack Opening (cm) 
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Figure 17. Crack Opening vs Time 


x 10> Relative Axial Displacement of Crack Surfaces 


Relative Axial Displacement (cm) 
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Figure 18. Relative Axial Displacement of Crack Surface 
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In order for the friction surface of a crack to 
dissipate energy, there must be both contact and relative 
motion between the surfaces of the crack. From the plots, it 
can be seen that contact is only experienced aa short 
time intervals, the relative axial motion during contact is 
small, and the amount of crack surface separation is also 
very small. Therefore, the damping ratio obtained through 
the logarithmic jeGEcHEnE equation is an averaged value of 
the nonlinear dissipation of energy experienced by the beam. 
Also, as the crack is placed further away from the 
cantilever wall, the amount of relative crack surface motion 
increases enexcay causing more. damping to occur. This 
explains the increase in damping ratio in the values of 


Table 8. 
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VII. CONCLUSIONS AND RECOMMENDATION 


The numerical es uihis of the embedded crack finite 
element model demonstrate the usefulness of the modeling 
technique in characterizing the dynamic response of a beam 
with interior through the thickness crack damage. 

Through the direct eanpeeisen of the different finite 
element model dynamic responses, it was shown that there was 
no concrete correlation between the amount of reduction in 
the flexural stiffness used in previous research and the 
response of an embedded crack finite element model. It was 
also found that the change in modal frequency experienced by 
the embedded crack model was quite small and was only 
comparable to a smaller reduced stiffness model vice the 30- 
50% reduction used in previous research. Furthermore, mode 
shape curvature was a more sensitive indicator of damage and 
a better locator of damage than modal frequency or mode 
shape displacement. The dizece Sopeoack in modeling an 
embedded crack was more accurate in capturing the true 
changes in dynamic response for the different modes of 
vibration. 

From the sandwich composite model analysis, it was 


shown that boundary conditions have a less significant 
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effect on the ability to detect damage from modal response 
changes. It was also shown that crack size was not directly 
proportional to the amount of change in modal response for 
all modes of vibration. However, increasing crack size does 
increase the change in the modal response. Crack location 
affects the ability to detect damage in certain modes of 
vibration. Therefore, evaluating a combination of different 
modes increases the chances of detecting damage from modal 
response. 

The nonlinear transient analysis furthered the study of 
damage detection by capturing the nonlinear crack opening, 
crack closure, and relative crack surface sliding motion 
during the transient response of a cantilever beam with an 
embedded crack. The friction environment between the 
surfaces of an embedded crack dissipates energy proportional 
to the amount of damping experienced by the beam during the 
transient response. The change in the amount of damping 
present can provide a ene. Oi seeee ns crack damage and 
Sonneouing the extent of crack damage present within the 
Structure. 

Recommendations for further study are as follows: 
investigate higher modes of vibration using the transient 
nonlinear analysis and frictional dissipation of energy, 


model multiple embedded cracks of various sizes’ and 
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locations, conduct experimentation with both homogenous and 
composite beams with physical damage and compare to results 


generate by numerical methods. 


45 








THIS PAGE INTENTIONALLY LEFT BLANK 


46 














LIST OF REFERENCES 


Adams, R.D., Cawley, P., Pye, C.d., and Stone B.J., 1978, “A 
Vibration Technique for Non-Destructively Assessing the 
Intergrity of Structures,” Journal of Mechanical Engineering 
Sciences, 20, 93-100. 


Kwon, Y. W. and Bang, H.C., 1997, The Finite Element Method 
Using MATLAB. CRC Press, Boca Raton, FL. | 


Pandey, A.K. and Biswas, M., 1994, “Damage Detection in 
Structures Using Changes in Flexibility,” Journal of Sound 
and Vibration, 169(1), 3-1/7. 


Petro, Samer H., Chen, Shen-En, GangaRao, Hota V.S., and 
Venkatappa, Suhas, 1997, “Damage Detection Using Vibration 
Measurements,” Proc. of the 15% International Modal 
Analysis Conference, Society of Experimental Mechanics, Il, 
dip Bes fo eI 


Rehm,G., 1987, “Nondestructive Tests for Early Damage 
Detection,” Proc. US-European Workshop on Bridge Evaluation, 
Repair, and Rehabilitation, St. Remy-Les-Chevreuse. 


Spyrakos, C., Chen, H.L., Stephens, J., and Govidara V., 
1990, “Evaluating Structural Deterioration Using Dynamic 
Response Characterization,” Proc. Intelligent Structures, 
Elsevier Applied Science, 137-154. 


Yuen, M.M.F., 1985, “A Numerical Study of Eigenparameters of 
a Damaged Cantilever,” Journal of Sound and Vibration, 103, 
301-310. 


4°7 














THIS PAGE INTENTIONALLY LEFT BLANK 


48 

















BIBLIOGRAPHY 


Bernasconi, Dr. O. and Ewins, Prof. D.J., “Application of 
Strain Modal Testing to Real Structures,” Proc. of the 7 
International Modal Analysis Conference, Society of 
Experimental Mechanics, 1453-1464. 


Ceravola, R. and De Stefano, A., 1995, “Damage Location in 
Structures Through a Connectivity Use of FEM Modal 
Analyses,” Modal Analysis: The International Journal of 
Analytical and Experimental Modal Analysis, Politecnico di 
Torino, 10, 178-186. : | 


Crespo, C., Ruotolo, C., and Surace, C., 13996, “Non-Linear 
Modeling of a Cracked Beam,” Proceedings of the 14™ 
International Modal Analysis Conference, Society of 
Experimental Mechanics, 1017-1022. 


Debao’, L., Hongcheng, Z., and Bo, Wang, 1989, “The 
Principles and Techniques of Experimental Strain Modal 
Analysis,” Proc. of the 7 International Modal Analysis 
Conference, Society of Experimental Mechanics, 1285-1289. 


Dong, C., Zhang, P.Q., Feng, W.Q., and Hang, T.C., 1994, 
“The Sensitivity Study of the Modal Parameters of a Cracked 
Beam,” Proc. of the 7 International Modal Analysis 
Conference, Society of Experimental Mechanics, 98-104. 


Gudmundson,P., 1982, “Eigenfregquency Changes of Structures 
Due to Cracks,” Journal of the Mechanics and Physics of 
Solids, 30(5), 339-353. 


Kam, T.Y. and Lee, T.Y., 1992, “Detection of Cracks in 
Structures Using Modal. Test Data,” Engineering Fracture 
Mechanics, 42(2), 381-387. 


49 














Lindner, D.K., and Kirby, G., 1994, “Location and Estimation 
of Damage in a Beam Using Identification Algorithms,” 
Proceedings of 35°" ATAA/ASME/ASCE/AHS/ASC Structures, 
Structural Dynamics and Materials Conference, 192-198, AIAA- 
94-1755-CP. 


Paolozzi, A. and Peroni, I., 1990, “Detection of Debonding 
Damage in a Composite Plate Through Natural Frequency 
Variations,” Journal of Reinforced Plastics and Composites, 
9; 369-309. 


Rizos, P.F., Aspragathos, N. and Dimarogonas, A.D., 1990,. 
“Identification of Crack Location. and Magnitude in a 
Cantilever Beam from the Vibrational Modes,” Journal of 
Sound and Vibration, Academic Press Limited, 138(3), 381- 
393% : 


Tsang, W.F., 1990, “Use of Dynamic Strain Measurements for 
the Modeling of Structures,” Proceedings of the eg 
International Modal Analysis Conference, Society of 
Experimental Mechanics, 1246-1251. 


50 








INITIAL DISTRIBUTION LIST 


Defense Technical Intormacton CONCSl™ sasieesens asinine 
8725 John J. Kingman Rd., STE 0944 
Ft. Belvoir, Virginia 22060-6218 


|B Y6Ye Wk 5am. o voy, quea Wwe of afb a' apere eon ne eR rarer nn menos Torte rrr eteri pr rrre ere erent rir sess 
Naval Postgraduate School | 
411 Dyer Rd. 

Monterey, California 93943 


Professor YOUNG W. Kwon, Code ME/ KW nncssssssssssssusssssssesecceseeeee wae 
Naval Postgraduate School 
Monterey, California 93943 


LT Stephen A. LipPSeSy, USN... ecsescsscsscscresscssasessnsnsesensossnenssosssesssnssssnsaessssess 
4 Indigo Court ; 
Savannah, Georgia 31406 


Naval Engineering Curricular Office, Code 34 vane - 


Naval Postgraduate School 
Monterey, CA 93943-5000 


ay § 


